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The pseudo-spin symmetry (PSS) is investigated in the density-dependent relativistic 
Hartree-Fock theory by taking the doubly magic nucleus ^■^^Sn as a representative. It is 
found that the Fock terms bring significant contributions to the pseudo-spin orbital poten- 
tials (PSOP) and make it comparable to the pseudo-centrifugal barrier (PCB). However, 
these Fock terms in the PSOP are counteracted by other exchange terms due to the non- 
locality of the exchange potentials. The pseudo-spin orbital splitting indicates that the PSS is 
preserved well for the partner states (i^3si/2, ^^2^3/2) of ^'^^Sn in the relativistic Hartree-Fock 
theory. 
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The pseudo-spin symmetry (PSS), quasi-degeneracy of two single particle states with the quan- 
tum numbers {n, l,j = l + 1/2) and (n = n — _L, Z_= l + 2,j=j = l + 3/2), was firstly discovered 

Jl . During the 
and the 

conservation conditions P, 0, liol . Illl . Il^ within the framework of relativistic mean field (RMF) 
theory. In Ref. 1^, the PSS was interpreted as the relativistic symmetry in the Dirac equation, 
which arises from the cancellation between an attractive scalar potential S5 and a repulsive vector 
potential Sq, i.e., Eg -|- Sq— 0, and the pseudo-orbit number / is nothing but the orbital angular 
momentum of the lower component of the Dirac wave function. However, in this limit there is no 
bound states in the mean field, at variance with reality. More realistic conditions were discussed 
for the conservation of the PSS in Refs. 



more than thirty years ago in spherical nuclei [l|, y| and later in deformed nuclei 
past ten years, many efforts have been made to investigate the origin of PSS |5 




let 

Q 



In the framework of the RMF theory 1^ 
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i.e., the relativistic Hartree 
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approach with a no-sea approximation, the nucleons interact via the exchange of mesons and 
photons. For the description of nuclear structme, the relevance of relativity is not in the need of 
relativistic kinematics but it lies in a covariant formulation which maintains the distinction between 
scalar and vector fields (more precisely, the zeroth component of the Lorentz four- vector field). The 
representations with large scalar and vector fields in nuclei (of the order of a few hundred MeV) 
provide more efficient descriptions of nuclear systems than non-relativistic approaches, for example 



the origin of the nuclear spin-orbit potential 



2l| and that of the PSS 



Although there exist some attempts to include the exchange terms in the relativistic description 



of nuclear matter and finite nuclei 



22 



23, l2J| , the relativistic Hartree-Fock (RHF) method was still 
not comparable with the RMF theory in the quantitative description of nuclear systems. Recently, 
it was shown that the density-dependent relativistic Hartree-Fock (DDRHF) theory [2^ gives a 
successful quantitative description of nuclear matter and finite nuclei at the same level as the RMF 
[2^. Compared with RMF, the relevance of the relativity is still kept well in DDRHF although 
the covariant formulation becomes much more complicated due to the exchange terms. Since the 
DDRHF describes quite well nuclear systems, it is worthwhile to investigate the role of exchange 
terms in the PSS, especially the influence of the non-locality on the conservation of the PSS. In this 
work, we study the role of the Fock terms on the conservation of PSS in the covariant relativistic 
Hartree-Fock theory. The numerical study is done with the effective Lagrangian PKOl [2^ in the 
case of DDRHF, and the results are compared with those obtained with the RMF model PKDD 

Q 

For spherical nuclei, the Dirac spinor can be written as. 



where xi{Ta) is the isospinor, Ga and Fa correspond to the radial parts of upper and lower compo- 

2 

nents, respectively. 'S^j^^a ^® spinor spherical harmonics and ^■^''ma(^) ~ ~^ ' ^'^jama^^) with 
= 2ja — la- In the spherical nuclei, the total angular momentum jaj its projection on the z axis 
nia and k = —(3{a- -L + l) are conserved. The eigenvalues of k are Ka = ±(ia + \) {— for ja = la + \ 
and -|- for ja = la — \)- In the following, the sub-index a will be omitted in the notations. 

Within the DDRHF, the radial Dirac equations, i.e., the relativistic Hartree-Fock equations for 



spherical nuclei, are expressed as the coupled differential-integral equations 



EG{r) 



d K 
dr r 



F{r) 



+ [S5(r) + So(r)]G(r)+y(r), 

EF{r) =+\^ + -] G(r) 
[dr r J 

- [2M + Ss(r) - So(r)] F(r) + X{r), 



(2a) 



(2b) 



where the scalar potential S5 and the time component of the vector potential Sq contain the 
contributions from the direct terms and the rearrangement term due to the density-dependence of 
the meson-nucleon couplings. The X and Y functions represent the results of the non-local Fock 
potentials acting on F and G, respectively {2^. By introducing the functions Xq, Xp, Yq and Yp 



as in Ref. 



_ G{r)X{r) nr)X{r) 

+ ^^^^ + + ^^"^ 

=XG{r)G{r)+Xp{r)F{r), 
=YG{r)G{r)+Yp{r)F{r), 



Y(r) 



(3a) 



(3b) 



the coupled differential- integral equations can be transformed into the equivalent local ones, 

F(r) - [A(r) - E] G{r) = 0, (4a) 
G(r) + [V{r) - E] F{r) = 0, (4b) 



Yp{r) 

dr r 



— + r 
ar r 



where A = + Yg, V = + Xp, and 



=Sn - Sc - 2M. 



(5a) 



27|. In 



The coupled equations (jlj can be solved by using the same numerical method as in RMF 
Eqs. (0J), the functions Xq, Xp, Yq and Yp might be taken as the effective potentials from the 
exchange (Fock) terms in the DDRHF. Eqs. @ must be solved iteratively since the potentials 
depend on the solution (G,F). 

^From the radial Dirac equations (0J), the Schrodinger-type equation for the lower component 
F is obtained as. 



-^F + Vi^F+ (VpcB + Vpso + V2)F 
dr'^ dr 

= _ (yi) _ E) (A^ - E) F, 



(6) 
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with 



Vi = {Xg - Yf) 



1 dA 



V2 = Yf 



1 (fA 
A-E dr 

D 



XgYf - ^Yf 
dr 



Vpso = 



+ Yg -E)+Xf{A-E), 

1 dA , 

{Xg + Yf) 



A-E dr 



k{1 — k) 



(7a) 

(7b) 
(7c) 

(7d) 



where VpcB and Vpso correspond to the pseudo-centrifugal barrier (PCB) and pseudo-spin orbital 
potential (PSOP), respectively. In the limit of Vpso = 0, the pseudo-spin becomes a good symmetry 
and the PSS can be labeled by the pseudo radial number n = n — 1, pseudo-orbit 1 = 1', and 
pseudo-spin s = s = ^, with the total angular momentum j = l±s for the two partner states. For 
instance, the partners are (nsi/2, {n — 1)^3/2) for / = 1, {np3/2, {n — l)/5/2) foi^ ^ = 2, etc. Notice 
that Xg and Yf in the Vpso are new contributions from the Fock terms compared with RMF. The 
potential V2 entirely originates from the Fock contributions. The direct (Hartree) and exchange 
(Fock) contributions of the PSOP and Vi can be separated as, 

K 1 dA^ 



D 

PSO 



r A — E dr 



D 



1 dYG 
A-E dr 

1 dA^ 

A-E dr ' 



Vf = {Xg - Yf) 



{Xg + Yf) 



1 dYG 



(8a) 
(8b) 

(8c) 
(8d) 



A-E dr ' 

while A comes from both the Hartree and Fock contributions. 

To have better a understanding of the pseudo-spin orbital splitting, especially the effects of 
Fock terms, it will be more transparent to rewrite Eq. (jEJ as, 



1 



d^ 



-E dr'^ 



F + 



1 



V^-E 



VpcB + + 



+ A^F = EE 



F 



(9) 



where the operators and are 



V 



D 



(10a) 
(10b) 



Let us discuss the PSS by taking the doubly magic nucleus ^^^Sn as a representative. In Fig. 
Uare shown the neutron single-particle energies calculated by the DDRHF with PKOl, compared 



with those given by the RMF with PKDD and the experimental data |23|- Among the neutron 
single-particle states shown in Fig. U there are four pseudo-spin partners Ip, Id, If and 2p, which 
correspond to the pairs {i^2si/2,i^ld^/2) , (i^2p3/2, 1^1/5/2) , (1^2^5/2,1^157/2) and {iy3si/2,i^2d-i/2) , 
respectively. From Fig. one can find a good PSS in the partners J^3si/2 and i^2d^/2 both in the 
DDRHF and RMF. 




PK01 Exp. PKDD 



FIG. 1: Neutron single particle energies of ^'^^Sn calculated by the DDRHF with PKOl and the RMF with 
PKDD. Experimental data are taken from Ref. 



^From the single-particle energies, the pseudo-spin orbital splitting is estimated as Aii^pso = 
{Ej-^i^^l^ — Ej-^jj^^i^) / {21 + 1) for the pseudo-spin partners Ip, Id, 1/ and 2p. The results are 
shown in Fig. [21 as a function of the average binding energy -Epso = (-^[j=[_i/2 + ^[j=/+i/2)/^' 
For both DDRHF (filled symbols) and RMF (open symbols) results, the pseudo-spin splitting 
between i^3si/2 and v2d^i2 is more than ten times smaller than that between z^2si/2 and 7^1^3/2. As 
shown in Fig. Uand Fig. |21 there exist some differences in the single-particle energies (especially 
the deeply bound states) between the DDRHF and RMF, whereas the monotonous decreasing 
behavior of AE'pso with the decreasing binding energies is observed in the both models. As a 
reference, the spin-orbit splitting AE'so = (^/j=z-i/2 ~ -^/j=/+i/2)/(2^ + 1) versus the average 
binding energy £^so = (^Zj+/-i/2 + ^/j=/+i/2)/2 is also shown in Fig. |2lfor the spin-orbit partners 
Ip, Id, 1/, and I5, and 2p, 2d. A clear difference between the pseudo-spin orbital and spin-orbit 
splittings can be seen in their energy-dependence, i.e., a strong energy-dependence is found for 
the pseudo-spin orbital splitting, while the spin-orbit splitting shows a weak energy-dependence. 
This difference can be understood from the comparison between the expressions for the PSOP (f7c|l 
and the corresponding spin-orbit potential. The energy E and the potential A in the denominator 
A — ii^ of Eq. (|7c() are comparable so that it brings a strong energy-dependence for the PSOP. On 
the other hand, the corresponding denominator of the spin-orbit potential in RHF is V — E ^, 
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which gives much weaker energy-dependence because the single-particle energy E is much smaller 
than the potential V: the value of is a few tens MeV while V is several hundred MeV. 

1.5 
1.2 
0.9 

> 

<1> 0.6 

s 

UJ 0.3 
0.0 
-0.3 

-50 -40 -30 -20 -10 

E (MeV) 

FIG. 2: The pseudo-spin orbital splitting Ai?pso — (-^/j=;%i/2 ~ ^ij=i+i/2) li"^^ + ^) versus the average 
binding energy -Bpso ~ (^/j=[-i/9 + -^rj=f+i/9)/2 for the neutron states in ^^^Sn. The pseudo-spin partners 
Ip, Id, If and 2p correspond to (i^2si/2, 1^1^3/2) , (j^2p3/2, 1^1/5/2) , (2^2^5/2,1^1.97/2) arid (1/351/2,2^2^3/2) 
states, respectively. The spin-orbit splitting A_Eso = (-£'/j=/-i/2 ~ £'j"=/+i/2)/(2^ + 1) are also given for 

(l/lp3/2, 1^1^1/2), (i/ld5/2,!^lc^3/2), (i^l/7/2, t^l/5/2), (j^l59/2 , ^^137/2) and (1^2^3/2,2^2^1/2), (2/2^5/2,2^2^3/2) 

pairs as a function of Ego = {Eij^i_i/2 +-£'y=(-Hi/2)/2- The results are obtained by the DDRHF with PKOl 
(filled symbols) and the RMF with PKDD (open symbols), respectively. 

Let us discuss the effects of Fock terms, especially the non-locality effect on the conservation of 
the PSS by considering the pseudo-spin partners Ip and 2p. The left panel of Fig. [^Jshows the radial 
wave function (upper and lower components G, F) for the pseudo-spin partners (z/2sx/2, 1^1(^3/2) 
and (1/851/2,2^2^3/2). The lower components of the partner states z/2si/2 (^^351/2) and u\d^j2 
(1/2^3/2) are very close each other in both the shape and magnitude whereas the upper components 
are quite different. Comparisons of the upper and lower components between the partner states 
indicate that the approximate PSS exists in the lower components for both partners Ip and 2p 
although the upper component of the total Dirac wave function dominate. The functions X{r) 
(y(r)) introduced in Eq.Q and representing the action of the Fock potentials on F (G) are also 
given in the right panel of Fig. I^for the pseudo-spin partners Ip and 2p. It is interesting to find 
that the radial dependence of the functions X and Y are very similar to those of the lower and 
upper components, respectively. 

In Eq. Q, the Fock-related terms X and Y are divided into two products of effective potentials 
times the Dirac wave functions to obtain the equivalent local Dirac equations (@J. Fig. |4l shows 
these effective potentials {Xq, Yq), and {Xp, Yp) for the pseudo-spin partners [u2si/2, and 
(z/3,Si/2, 1/2^3/2). The local peaks and dips in the figures are due to the localization of the exchange 
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FIG. 3: The radial wave functions G and F (left panel), and the nonlocal terms X and Y (right panel) in 
Eq. (|2Jl given by the DDRHF with PKOl for the pseudo-spin partners (i/2si/2, ^^1^3/2) ^ind (j^3si/2, ^^2^3/2) 



in i32Sn. 

terms in Eq. reflecting the nodes of the upper components in the denominator. It is seen in 
Fig. EJthat the radial dependence of Yp is almost identical to that of Xg since the shapes of {X, Y) 
and (G, F) are very close each other. Because the peaks and dips in Fig. EJare due to the nodes of 
upper components, their contributions to the single-particle energies will be smoothed out by the 
corresponding nodes. Thus, one can ignore the non-locality represented by the peaks or dips for 
the discussion on the PSS. Besides the peaks, one can find that in the inner part of the nucleus the 
effects of the non-locality are significant on Xq and Yp and give a strong state dependence with 
different signs, while it shows fairly weak effects on Yq. The Xp is just due to the nodes of the 
upper component. 
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FIG. 4: The effective potentials Xq, Yq (left panel), and Xp, Yp (right panel) due to the exchange (Fock) 
terms in Eq. for the pseudo-spin partners (z^2si/2, J^ld3/2) and (1^351/2, 1^2^3/2) in ^■^^Sn. See the text 
for detail. 



In RMF it has been proved that the PSS is well obeyed if the PCB is much stronger than the 
PSOP la]. Fig. Elshows the PCB and PSOP multiplied by the factor F^/ (y^ - E) for the pseudo- 
spin partners (2^251/2)^^1^^3/2) and (1^351/2,^^2^3/2) in "'^'^^Sn. Due to the denominator A — in Eq. 
(|7c|) . there exists singular points at r ~ 6 fm for the partner (i/2si/2, 2^1(^3/2) , and at r ~ 7.5 fm 
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for (z^3si/2i ^^2^3/2) • The other local peaks in the PSOP are due to the nodes of upper component 
G. For the s states (/ = 0), the PCB is much stronger than the PSOP since the contributions for 
the PSOP around the nodes or the singular points are more or less mutually cancelling. On the 
other hand, for the d states, the PSOP are comparable to the PCB even after taking account of 
the cancellation around the nodes or singular points. Comparing the shapes of the PSOP in Fig. 
ISl and Xq , Yp in Fig. |4l one can find that the Fock terms present significant contributions to the 
PSOP, especially for the d states in the inner part of the nucleus. It is also seen that the Fock 
terms in Fig. ISJhave substantial contributions to the PSOP in Fig. |5| 




FIG. 5: The PCB and PSOP multiplied by the factor F^/ (l/^ - E) for the pseudo-spin partners 
{y2sif2Tv\dy,/2) and {v'isi/2iv2d'^/'2) ^'^^Sn. The PCB contributions are drawn by the dashed lines, 
while the PSOP are shown by the solid lines. The DDRHF with PKOl is used for the calculations. See the 
text for details. 



^From Eq. (jUJ, one can estimate the contributions of the potentials I^cBj a-^d to the 
single-particle energy E. For example, the PCB contribution can be expressed as, 



J [°° ^PCB , 

F^drJo VD-E 



(11) 



The results calculated by the DDRHF with PKOl and the RMF with PKDD are shown in Table 
Q] for the pseudo-spin partners Ip and 2p. For both the DDRHF and RMF, the terms F" , /SP 
and in Eq. Q show substantial differences between the partner states (z^2si/2, '^1^3/2) and 
(z^3si/25 '^2^3/2) whereas the differences in the PCB and the exchange terms are negligible. The 
differences in F" and reflect those of the lower components in the two pseudo-spin partners as 
shown in the left panel of Fig. |31 

The large differences between the partner states can be seen in the F" ^ and contributions. 
However, these three terms cancel largely one another and the PSS is preserved to a good degree, 
especially in the partners (i^3si/25 both DDRHF and DDRMF. The Fock terms 

become small although each term on the right hand side of Eq. ()10b|) shows appreciable differences 
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TABLE I; The single-particle energies E and the eontributions from different terms in the left hand side of 
Eq. © given by the DDRHF with PKOl, in comparison with those by the RMF with PKDD. AU units are 
in MeV. 



Model 


Orbit 


E 


F" 


^PCB 






DDRHF 
PKOl 


I/2si/2 

i^ldg/a 


-31.41 
-34.90 


18.11 -75.35 
14.87 -79.01 


9.30 
9.54 


-2.99 
0.44 


19.51 
19.26 


t^3si/2 


-8.33 
-8.66 


34.25 -72.00 
31.93 -73.96 


11.11 
11.32 


0.09 
3.89 


18.22 
18.17 


DDRMI 
PKDD 


I/2.S1/2 
Vld3/2 


-34.81 
-38.87 


21.86 -64.65 
18.17 -68.08 


11.04 
11.41 


-3.07 
-0.37 




i^3si/2 

1/24/2 


-8.15 
-8.44 


40.13 -61.97 
37.65 -63.75 


13.02 
13.36 


0.67 
4.30 
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FIG. 6: The fimctions V^p^o^V {V^ - ^) and VfFF'/ {V^ - E) given by the exchange (Fock) terms of 
the DDRHF with PKOl for the pseudo-spin partners {v2sii2TV^d^/'2) and (1^351/2,^^2^3/2) of ^■^^Sn. The 
singular points for z^2.Sj/2 and vld'^i2 at r ~ 6 fm. and for ;/3si/2 and vldj^j^ at r ~ 7.5 fm are due to the 
denominator (A — E) in the PSOP (see Eq. lf7c|l ') while the other local peaks are due to the nodes of the 
upper component G (see Eq. 

between the partner states: for example, large for d-states. The contributions from the 

exchange potentials and in Eq. Q are shown in Fig. El For the s states, the exchange 

terms and give small contributions because of their changing signs. On the other hand, 

for the d states there are significant cancellations between and ■, especially in the inner 

part of the nucleus. Although the Fock terms bring substantial contributions to the PSOP, these 
contributions are cancelled by the other exchange term Vy ^ which stems mainly from the non- 
locality (the state-dependence) of the exchange potentials. Thus, the PSS still remains preserved 
even after the inclusion of Fock terms due to these large cancellations as can be seen in Table Q] 
Let us now discuss the reason why the cancellations among the exchange terms occur. From 
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the similar radial dependence between the non-local terms X iY) and Dirac wave functions F (G) 
in Fig. El we might be able to validate the following relations, 

X{r) ~Xo(r)F(r), Y{r) c^Yo{r)G{r), (12) 

where Xq and Yq are supposed to be state-independent potentials due to the Fock terms. Using Eq. 
H12|) . the non-local RHF equations can be reduced to local ones similar to the RMF equations 
in which the terms Xq and Yp do not appear in the PSOP. Thus, the realization of the PSS will 
be similar to the RMF case. Therefore, the cancellation of the Fock contributions in Table Q] is 
not occasional but it is because of the similar radial dependence between the Fock-related terms 
{X,Y) and the wave functions (F, G). 

In summary, the PSS in the DDRHF theory was investigated in the doubly magic nucleus ^'^^Sn. 
The PSOP was derived by transforming the coupled radial Dirac equations into the Schrodinger 
type equation of the lower component properly taking account of the non-local Fock terms. The 
analyses of the single particle spectrum and the pseudo-spin orbital splitting indicate that the 
PSS is preserved as a good symmetry for the pseudo-spin partner (z/3si/2) ^^2^3/2) of -"^^^Sn in the 
DDRHF on the same level as RMF, although the Fock terms bring substantial contributions to 
the PSOP. These contributions to the pseudo-spin orbital splitting. However, are cancelled by the 
other terms due to the non- locality of the exchange potentials. The physical mechanism of these 
cancellations was discussed in relation to the similarity between the exchange potentials and the 
Dirac wave functions. 
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